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Extremely high precision of near-future radio/optical interferometric observatories 
like SKA, Gaia, SIM and the unparalleled sensitivity of LIGO/LISA gravitational- 
wave detectors demands more deep theoretical treatment of relativistic effects in 

>, 

the propagation of electromagnetic signals through variable gravitational fields of 

VO . 

the solar system, oscillating and precessing neutron stars, coalescing binary sys- 



o 
o 

cr 



terns, exploding supernova, and colliding galaxies. Especially important for future 
gravitational-wave observatories is the problem of propagation of light rays in the 
field of multipolar gravitational waves emitted by a localized source of gravitational 
radiation. Present paper suggests physically-adequate and consistent solution of this 
problem in the first post-Minkowskian approximation of General Relativity which 

x. 

accounts for all time-dependent multipole moments of an isolated astronomical sys- 
tem. We derive equations of propagation of electromagnetic wave in the retarded 
gravitational field of the localized source emitting gravitational waves of arbitrary 
multipolarity and integrate them analytically in closed form. We also prove that 
the leading terms in the observable relativistic effects (time delay, deflection angle, 
and rotation of the plane of polarization of light) depend on the instantaneous value 
of the multipole moments of the isolated system and its time derivatives taken at 
the retarded instant of time but not on their integrated values. The influence of 
the multipolar gravitational field of the isolated system on the light propagation is 
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examined for a general case when light propagates not only through the wave zone of 
the system but also through its intermediate and near zones. The gauge freedom of 
our formalism is carefully studied and all gauge-dependent terms are singled out and 
separated from observable quantities. We also present a thorough-going analytical 
treatment of time delay, light-ray bending and polarization of light in the case of 
large impact parameter corresponding to the approximation of a plane gravitational 
wave of arbitrary multipolar ity This exploration essentially extends previous results 
regarding propagation of light rays in the quadrupolar field of a plane monochromatic 
gravitational wave. Explicit expressions for time delay, deflection angle and rotation 
of the plane of polarization of light are obtained in terms of the transverse-traceless 
(TT) part of the space-space components of the metric tensor. We also discuss the 
relevance of the developed formalism for detection of relativistic effects of multipolar 
gravitational fields by existing and near-future astronomical techniques as well as 
the gravitational wave detectors. 
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1. INTRODUCTION 

1.1. Statement of the Problem and Relation to Previous Work. 

We consider propagation of electromagnetic signals through the time-dependent grav- 
itational field of an isolated gravitating system emitting gravitational waves such as, for 
example, a binary system or oscillating neutron star. Gravitational field is assumed to be 
weak everywhere along the light ray and a linear (post-Minkowskian) approximation of Gen- 
eral Relativity is used. We do not restrict our consideration to a specific system, but rather, 
following [1] , develop a formalism which can be applied to an arbitrary self-gravitating source 
that is completely characterized by its time-dependent multipole moments of arbitrary or- 
der. We assume the wavelength of light to be much smaller than the wavelength of the 
gravitational waves which allows us to work in the geometrical optics approximation. The 
relative distances between and positions of the source of light, the isolated system, and the 
observer are not restricted as well, which makes our mathematical formalism quite general 
and applicable for most practical situations. With the assumptions made, we show that 
the equations of light propagation, which in geometrical optics are the equations of null 
geodesies, can be analytically integrated giving linearized metric perturbations to the un- 
perturbed flat space-time trajectory. We also use the formalism developed in [2] and [3] to 
consider the propagation of the polarization of light along the light ray and, in particular, 
the effect of gravity-induced rotation of the polarization plane of an electromagnetic wave 
(the Skrotskii effect). 

This work is essentially a generalization of [1] for the case when the isolated system pos- 
sesses multipole moments of arbitrary order and type. In [1] the same problem is considered 
in spin-dipole and mass-quadrupole approximation only. We also generalize the results of 
[3] for Skrotskii effect on the case of source possessing higher order multipole moments. The 
results of [3] for isolated system emitting gravitational waves were obtained again only in 
spin-dipole mass-quadrupole approximation. Formalism of the present paper allows us to ex- 
tend results of exploration of propagation of electromagnetic waves in random field of plane 
gravitational waves [4, 5] to more general class of an isolated systems emitting non-planar 
gravitational waves anisotropically. Additional mathematical extention of the present work 
with more thorough historical review isgiven in the next publication [6]. 
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1.2. Notations and Conventions 

Notations in this paper coincide with those of [6]. Metric tensor on the space-time 
manifold is denoted by g a p and its perturbation h a p = g a p — r] a p. The determinant of 
the metric tensor is negative and is denoted as g — det||<7 a( g||. The four-dimensional fully 
antisymmetric Levi-Civita symbol e a ^& is defined in accordance with the convention £0123 = 
+ 1. In the present paper we use a geometrodynamic system of units [7, 8] such that c = 
G — 1 where c is the fundamental speed and G is the universal gravitational constant. 
Space-time is assumed to be asymptotically-flat and covered by a single coordinate system 
(x°, x 1 , x 2 , x 3 ) = (t, x, y, z) , where t and (x, y, z) are time and space coordinates respectively 
This coordinate system is reduced at infinity to the Minkowskian coordinates. We shall also 
use the spherical coordinates (r, 9, <p) related to (x, y, z) by the standard transformation 

x = rsin#cos0, y = r sin # sin , z = rcos9. (1.1) 

Greek (spacetime) indices range from to 3, and Latin (space) indices run from 1 to 3. 
If not specifically stated the opposite, the indices are raised and lowered by means of the 
Minkowski metric r\ a p = diag(— 1, 1, 1, 1). Regarding this rule the following conventions for 
coordinates hold: x l = Xi and x° = —Xq. We also adopt notations, 5^ = diag(l, 1, 1), for 
the Kroneker symbol (a unit matrix), and, e^, for the fully antisymmetric 3-dimensional 
symbol of Levi-Civita with convention 6123 = +1. 

Repeated indices are summed over in accordance with the Einstein's rule [7, 8]. In the 
linearized approximation of general relativity used in this work there is no difference between 
spatial vectors and co-vectors or between upper and lower space indices. Therefore, for a 
dot product of two space vectors we have 

A l B % = AiBi = A 1 B 1 + A 2 B 2 + A 3 B 3 . (1.2) 

In what follows, we shall commonly use spatial multi-index notations for Cartesian three- 
dimensional tensors [9], that is 

%a, = 2oi...o,- (1-3) 

Tensor product of I identical spatial vectors k l will be denoted as a three-dimensional tensor 
having / indices 

kai...k at = k ai ai . (1-4) 
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Full symmetrization with respect to a group of spatial indices of a Cartesian tensor will be 
distinguished with round brackets put around the indices 

Q(oi...o,) = jj ^ Qg(i)...g(Q » 

where a is the set of all permutations of (1, 2, /) which makes Q ai ... ai fully symmetrical in 
a 1 ... a x . 

It is convenient to introduce a special notation for symmetric trace-free (STF) Cartesian 
tensors by making use of angular brackets around STF indices. The explicit expression of 
the STF part of a tensor Q ai ... ai is [9] 

[1/2] 

Q<ai...ai> = ^ ] a fc^(aia 2 ' ' ' ^a2k-ia2k^o, 2 k+i—ai)b 1 b 1 ...b k b k > (1-6) 
fc=0 

where [//2] is the integer part of 1/2, 

Sai...ai = Q(ai...aj) (1-7) 

and numerical coefficients 

g (2Z-2fc-l)!! 
fc (2jfe)!!(2Z-l)!! (Z-2A;)! ' 1 } 

We also assume that for any integer / > 

l\ = l(l- 1) . . .2 - 1 , 0! = 1, (1.9) 

and 

Z!! = Z(Z-2)(Z-4)...(2orl), OH = 1 . (1.10) 

One has, for example, 

1 1 1 

T <a bc> = 2~(a6c) — -^>abT( c jj) ~ ^bcT(ajj) — ^a C T(bjj) ■ (l-H) 

Cartesian tensors of the mass- type 1<a l > and spin-type multipoles <S<A;> entirely describing 
gravitational field outside of an isolated astronomical system are always STF objects that can 
be checked by inspection of the definition following from the multipolar decomposition of the 
metric tensor perturbation h a p [9]. For this reason, to avoid appearance of too complicated 
index notations we shall omit in the following text the angular brackets around indices of 
these (and only these) tensors, that is we adopt: = 1<a 1 > and iS^ ; = 5<a ; >- 
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We shall also use transverse (T) and transverse-traceless (TT) Cartesian tensors in our 
calculations [8, 9]. These objects are defined by making use of the operator of projection 
Pj k = 8jk — kjk onto the plane orthogonal to a unit vector kj. Thus, one has [9] 

Qa\...a x — Pa, 1 b 1 Pa 2 b 2 ---Pa l biQb 1 ...bi , (1-12) 
[1/2] 

Qai...ai = S ^ l ^ > kP(axa2 ' ' ' Pa 2k _ 1 a 2k Wa 2k+1 ...ai)b 1 b 1 ...b k b k ■, (1.13) 
fe=0 

where again [1/2] is the integer part of 1/2, 

W ai ... ai = QL ai) (1.14) 



and numerical coefficients 



■ (-!)*<(< -fc-l)H 
fc 4 fc fc!(Z-2ife)! ' 



;i.i5) 



For instance, 



= 2 (PaiPbjQij + PbiPajQij) — "-Pat {PjkQjk) ■ (1-16) 



Polynomial coefficients will be used in some of our equations and they are defined by [10] 

CKpi,...,^)^-^ r, (1.17) 

Pi! ...p n \ 

where I and pi are positive integers such that $^=1 ^* = ^ ^ e introduce a Heaviside unit 
step function H(p — g) such that on the set of whole numbers 

0, if p < q, 

H(p-q) = { ~ (1.18) 

1, if p > q. 

For any differentiable function / = f(t, x) one uses notations: / j0 = df/dt and f ti = df/dx l 
for its partial derivatives. In general, comma standing after the function denotes a partial 
derivative with respect to a corresponding coordinate: f, a = df(x)/dx a . Overdot denotes 
a total derivative with respect to time / = df /dt — df /dt + x l df /dx l usually taken in this 
paper along the light ray trajectory x(t) = (x l ). Sometimes the partial derivatives with 
respect to space coordinate x 1 will be also denoted as <9j = d/dx l , and the partial time 
derivative will be denoted as d t = d/dt. Covariant derivative with respect to the coordinate 
x a will be denoted as V a . 
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We shall use special notations for integrals with respect to time and for those taken along 
the light-ray trajectory. Specifically, the time integrals from a function F(t,x), where x is 
not specified, are denoted as 



t 



^«,.) S /F(r,-)*, F^^fp^r,^. (1.19) 

— oo — oo 

Time integrals from the function F(t, x) taken along the light ray so that spatial coordinate 
x is a function of time x = x(t), are denoted as 

t t 
F^(t,x)= J F(r,x(r))dr, F^(t,x) = J F^{r,x{r))dr . (1.20) 

— oc — oo 

Integrals in equations (1.19) represent functions of time and space coordinates. Integrals in 
equations (1.20) are defined on the light-ray trajectory and are functions of time only. 
Multiple time derivative from function F(t, x) is denoted by 

*»<*..) = ^f>, (1.21) 
so that its action on the time integrals eliminates integration in the sense that 

( ' >~ dtr+> ~ dtp*'' ' 1 ' 

and 

r„i, x dP +1 F^ l Ut,x) dP +2 F^ 2 Ut,x) 
F®(t, x) = , J - 1 = , J ' ; . 1.23 

Spatial vectors will be denoted by bold italic letters, for instance A 1 = A, k l = k, etc. 
The Euclidean dot product between two spatial vectors a and b is denoted with dot between 
them: a l bi = a ■ b. The Euclidean wedge (cross) product between two vectors is denoted 
with symbol x, that is €ijka^b k = (a x b)\ Other particular notations will be introduced as 
they appear in the text. 

2. METRIC TENSOR AND COORDINATE SYSTEMS 
2.1. Metric Tensor 



We assume that the gravitational field is weak everywhere along the light-ray path. This 
allows us to work in the linear (with respect to G) approximation of General Relativity which 



9 



is sometimes called the post-Minkowskian approximation [11]. Thus, the metric tensor can 
be expressed as the sum of the Minkowski metric and a small perturbation h a p <C 1 being 
proportional to the universal gravitational constant G: 



9 a/3 — Vaf3 + h af3 . 



(2.1) 



We consider propagation of light in the space-time around an isolated gravitating system. 
The most general expressions for metric perturbations in this case were given by Thorne [9] 
in a special coordinate system which belongs to the class of the harmonic coordinate systems. 
In what follows we shall call these special coordinates canonical harmonic coordinates. The 
expressions given by Thorne read 



Lean. 
"-00 



Lean. 
a 0i 



2M 



2 £- 

1=2 
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^^tpq^p 1 v q 
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/=2 
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(/ + !)! 



i S 3)pAi-2 (t ~ r ) 



(2.4) 
(2.5) 



Here .M and <Si are the total mass and angular momentum of the system, X\ l and Sa { are 
two independent sets of mass-type and spin-type multipole moments, and N l = x % jr is a 
unit vector directed from the origin of the coordinate system to the field point. Since the 
origin of the coordinate system has been chosen at the center of mass, the expansions (2.2) 
- (2.5) do not depend on the mass-type linear multipole moment = 0. 

The expressions for the metric perturbations in a different coordinate system can be 
obtained by means of an infinitesimal coordinate transformation 



x 



la 



x a -w a 



(2.6) 



from the canonical harmonic coordinates x a to coordinates x' a where w a are gauge functions. 
Then in the new coordinate system one has [8] 



h a p = Kp' + w a,(3 + tCfto • 



(2.7) 
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We note that if functions w a satisfy the homogeneous wave equation 



Uw a = 0, 



(2.8) 



the transformation (2.6) leaves us within the class of the harmonic coordinate systems, since 
the harmonic gauge conditions can be formulated as \Z\x a = [12]. 

The most general solution of the equation (2.8) is given in [9] (see also [13]) and reads 

'WAXt-r) 



°=E 

1=0 

oo 

1=0 

oo 

E 



,Al 



X Al (t-r) 



+ 



(2.9) 
(2.10) 



Mi 



yi Al _At-r) 



+ 



where Wa^ Xa 3^a ; _!, and 2 q A l _ 1 are four sets of STF multipole moments depending on 
the retarded time. 

We shall use the gauge freedom in choosing functions w a to simplify the problems of 
integration of the equations of light propagation and analysis of the observable effects. 
It turns out that two coordinate systems possess properties useful for solution of these 
problems. These are the harmonic and the Arnowitt-Deser-Misner (ADM) coordinates. 



2.2. The Harmonic Coordinate System 

The harmonic coordinate system, analogous to the Lorentz gauge in electrodynamics, is 
defined in linear approximation by the conditions [12] 

2h a ^ -h' a = (2.11) 

imposed on the metric perturbations, where h = h a a . The linearized Einstein field equations 
in harmonic coordinates reduce to wave equations for metric perturbations h a ^. 

The metric tensor perturbations in harmonic coordinates (2.2) - (2.5) possess a property 
that can be used to simplify the problem of analytical integration of the equations of light 
propagation. Specifically, the method of integration of equations of light geodesies, developed 
first in [15] for stationary fields, and extended in [1] to time-dependent fields, is utilized in 
the present work (Section 44.1). This method takes advantage of the fact that the metric 
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perturbations are functions of the retarded time. Not all coordinates in general relativity 
have this property. For instance, metric perturbations in the canonical ADM coordinate 
system, considered in the next section, contain terms that depend on instantaneous time, 
which does not allow us to use the abovementioned method of integration with the retarded 
time. 

2.3. The Arnowitt-Deser-Misner Coordinate System 

The ADM gauge conditions [14] in the linear approximation read 

2/ioi,j — hufi = , — hjj t i = . (2-12) 

For comparison, the harmonic gauge conditions (2.11) in the linear approximation have the 
form 

2^oi,j — Ki,o — hoofi , 2/iyj — hjj ti = — hooj . (2-13) 

The ADM coordinates have a property that test particles perturbed by gravitational 
waves stay at rest with respect to these coordinates, that is the ADM coordinate system is 
co-moving with test particles. This property can be used to simplify analysis of observable 
effects. For example, when one considers detection of a photon by an observer in the field 
of gravitational wave, then normally one would have to solve both the equations of motion 
of the photon and observer [12]. Using ADM coordinates allows us to exclude the problem 
of motion of the observer from the coordinate description of the gravitational-wave effect on 
the photon. 

The disadvantage of using the ADM coordinates for analysis of the effects in propagation 
of light through gravitational field of an isolated source is that, as it was mentioned above, 
the ADM metric perturbations normally contain terms that are functions of instantaneous 
time. This does not allow one to use mathematical advantages of the method of integration 
of the equations of light propagation (Section 44.1). 

Fortunately, the classes of ADM and harmonic coordinates overlap. In [1] an ADM- 
harmonic coordinate system was constructed in the case when isolated system possesses 
constant mass and angular momentum and time-dependent quadrupole moment. We gen- 
eralize this result to the case of isolated system possessing the whole multipole structure 
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(2.2)-(2.5). The gauge functions 



oo , y 



W 



W 




4 E 



^ (/ + !)! 



Z=2 



UbaS^ ^t - r) 



bring the metric perturbations h a p to the form 

2M 



h 



oo 



+ 



Ai-i 



(2.14) 
(2.15) 



(2.16) 



h 



o? 



^^ipqSpN q 



hij — 5ijh()Q + /l^- , 



(2.17) 
(2.18) 



UTT _ p can. 
"ij _ r ijkl ( ikl i 



(2.19) 



which satisfies both harmonic and ADM gauge conditions. Here the TT-projection dif- 
ferential operator Pijki, applied to symmetric tensors depending on both time and spatial 
coordinates, is given by 

P ijkl = (8 ik - A" 1 ^)^ - A^djdi) - i(c% - A^d^H - A^dkdi) , (2.20) 

and A and A -1 denote the Laplacian and the inverse Laplacian respectively. 

The ADM-harmonic coordinates combine the advantages of both coordinate systems and 
will be used in what follows for analysis of observable effects in propagation of electromag- 
netic signals through the gravitational field of the isolated system. 

3. EQUATIONS FOR PROPAGATION OF ELECTROMAGNETIC SIGNALS 



3.1. Geometrical Optics in Curved Space-Time 



In this section following [3] and [8] (see §22.5) we provide a derivation of the main laws 
of geometrical optics in curved space-time from the Maxwell's equations. 
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In absence of sources the Maxwell's equations for the electromagnetic field tensor F a/3 in 
curved space-time take the well known form [7, 8] 

V a F Pl + V p F ia + V 7 F Q/3 = , (3.1) 

V p F aP = , (3.2) 

where V a denotes a covariant derivative. The wave equation for F a p can be derived from 
(3.1) and (3.2) and in vacuum (R a /3 = 0, where R a p is the Ricci tensor) it assumes the form 

U g F aP + R a ^F^ = , (3.3) 

where D g = V Q V a and R a /3j8 is the Riemann curvature tensor of the space-time. 

The solution to the Maxwell's equations corresponding to a high frequency wave is given 
by [4, 5] 

F afS = Re{A afS exp(i(p)} , (3.4) 

where A af3 is the complex amplitude, which is a slowly varying function of position and time, 
and ip is the phase, rapidly changing with position and time. 

The criterion for applicability of the geometrical optics approximation in curved space- 
time is formulated [8] as follows [8]. Let 1Z be the characteristic radius of curvature of the 
space-time through which the electromagnetic wave propagates and C - the characteristic 
length, over which the amplitude, wavelength and polarization of the electromagnetic wave 
change significantly. Geometrical optics approximation can be applied whenever the wave- 
length of the electromagnetic radiation satisfies the following two conditions: 1)\<^H] and 
2)A <C C. Then a small parameter e = X/min{C, TZ} can be introduced, and the expansion 
of the electromagnetic field of a wave (3.4) in powers of e is assumed to have the form 

Fap = (a a p + eb a/3 + e 2 c a/3 + ...) exp ■ ( 3 - 5 ) 

For further discussion of this procedure see [16] and [19]. Substituting this expansion into 
(3.1), taking into account the definition of the electromagnetic wave vector l a = dip/dx a 
and rearranging the terms with the same powers of e lead to the chain of equations 



l a apj + a^ a + / 7 a a p — , 

V a Og 7 + V,3 a ja + V 7 a af 3 — — i(l a feg 7 + 1/3 b ya + / 7 b a p) , 



(3.6) 
(3.7) 
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where the effects of curvature appears in the next approximation and, hence, have been 
neglected. Another chain of equations is obtained by substituting (3.5) into (3.2): 

lp a aP = , (3.8) 
a a(B + Up b a/3 = . (3.9) 

The equation (3.8) shows that in the lowest order approximation (terms of the order l/e) 
the amplitude of the electromagnetic field tensor is perpendicular to the wave vector in the 
four dimensional sense. Contracting (3.6) with l a and accounting for (3.8) implies that l a is 
a null vector 

lj a = 0. (3.10) 

Taking a covariant derivative from the equality (3.10) and using the fact that V[ a lp] = 
(since l a = V a y?) we obtain that l a satisfies the geodesic equation 

fy p F = , (3.11) 

showing that the wave vector is parallel transported along itself. Equations (3.10) and 
(3.11) combined constitute an important result of the geometrical optics: in the lowest 
order approximation light rays are null geodesies. 

Substituting (3.4) into (3.3) and considering terms of the order l/e one can obtain the 
law of propagation for the amplitude of the electromagnetic tensor 

D\a a p + -&a a p = , (3.12) 

where D\ = l a V a and i? = (1/2) V^P is the expansion of the null congruence l a . 

3.2. Equations for Trajectory of a Photon 

It was shown in the previous section that in the lowest order approximation of geometrical 
optics light rays, or trajectories of photons, are represented by null geodesies of the space- 
time under consideration. In what follows we shall consider propagation of photons in 
a space-time around an isolated gravitating system. With the metric of the space-time 
defined by (2.7). We shall assume that the wavelength of the electromagnetic radiation is 
much smaller than the characteristic wavelength of the gravitational waves emitted by the 
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gravitating system, so that the conditions for validity of geometrical optics approximation 
are satisfied. 

Let us consider propagation of photons subject to the initial-boundary conditions 

/ \ dx(-oo) , ,„ „ „. 

x(to)=x , y dt J = k , (3.13) 

which specify the spacial velocity of photons at the past null infinity by a unit spacial vector 
k and the positions of the photons at some initial instant of time t . 

Taking into account that l a = dx a /d\, where A is an afline parameter along the photon's 
trajectory, we can rewrite the equations of geodesies (3.11) in terms of coordinates of a 
photon 

iv + ^ixix =0 ' < 3 - 14 > 

where 

r^ 7 = \g a " (drftf + - dtfpy) (3.15) 

are the Christoffel symbols. Combining the four equations (3.14) we obtain the equations 
for spacial position of a photon as a function of the coordinate time 

m = - (r^ - rV) . (3.16) 

Substituting (3.15) into (3.16), taking into account that x % = k l + 0(h) and keeping only 
linear in metric perturbations terms we can rewrite the equations (3.16) in the form 

^(t) =^ h oo,i - Vo - ^hoo,ok l - h ikfl k k - (h 0i , k - h ok ,i)k k - (3.17) 

hoo^k* - (h ikJ - hi kj ^j k k k j + Q/%, - h ok ^ k k k j k l . 

In section 4 we will show how these equations can be analytically integrated in the space-time 
with metric (2.7) defining gravitational field of an isolated gravitating system. 



3.3. Equations for Propagation of Polarization Properties 

In this section, following [2] and [3], we introduce the relativistic description of polarized 
electromagnetic radiation and give the equation (derived in [3]) for gravity- induced rotation 
of the polarization plane of an electromagnetic wave (the Skrotskii effect). 

For description of the polarization properties of electromagnetic radiation it is necessary 
to introduce local reference frames along the light rays. At each point of space-time we 
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introduce a complex null tetrad (l a ,n a ,m a ,m a ) associated with a congruence of light rays 
[20, 21]. All the vectors of the tetrad are lightlike. The vectors l a and n a are real; vectors 
m a and m a are complex and complex conjugate with respect to each other. The tetrad 
is normalized in such a way that the only nonvanishing products among the vectors are 
l a n a = — 1 and m a fn a = 1. 

It is obvious that such field of complex null tetrads is not uniquely determined by the 
congruence of null rays and the normalization conditions. The transformations 

l' a = Al a , (3.18) 
A' 1 (n a + Bm a + Bm a + BBl a ) , 
e~ ie (m a + Bl a ) , 
m' a = e ie (m Q + Bl a ) , 

where A and are real and B is a complex parameter, preserve the direction of the vector 
l a and the normalization conditions [21]. The transformations (3.18) form a four-parameter 
subgroup of the Lorentz group. In addition to the complex null tetrad we introduce at each 
point in space-time an orthonormal reference tetrad e"^ defined as follows. Suppose at each 
point of space-time there is an observer moving with four-velocity u a . Let two of the vectors 
of the observer's local reference frame be 

e Q (0) = u a , e a {3) = (-I/)" 1 [r + (l p u p )u a ] . (3.19) 

With such orientation of the reference frame the observer will see the electromagnetic wave 
propagating in the +z direction. Spacelike unit vectors e"^ and e a ^ are orthogonal to each 
other as well as to both e"^ and e a ^ and thus are specified up to a spacial rotation. 

The relationship between the vectors of the complex null tetrad and the frame e"^ is 
given by 

r = -M( e ; ) + e y , (3.20) 
i 

2 

1 

72 

"<■ —('":• "V) • (3-23) 



a = -^){e a (0) -e a {3) ) , (3.21) 



nr =^=(e a {1) +ie a {2) ) , (3.22) 



V2 
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Vectors e a (i) and e a ( 2 ) and corresponding to them m a and m a play an important role in 
description of polarization properties of electromagnetic radiation since they form a basis in 
the polarization plane. 

The tensor of the electromagnetic field can be represented as 

F a(3 = Re{F aP ) , (3.24) 

where T a p is the complex field. The complex field can be expressed as 

Fap = <&l[amp] + ^l\0fn a] (3.25) 

where <£> and \1/ are complex scalar functions. 

The electric and magnetic fields in the rest frame of an observer moving with four- velocity 
u a are defined as 

E a = F a % , H a = (-l/2)e a/37<5 F 75 M /3 . (3.26) 

We also define the complex electric field as 

S a = T^up . (3.27) 

Polarization properties of light are completely characterized by the polarization tensor [7] 
(coherency matrix in [22]) 

J a p = (£ a £p) , (3.28) 

where the angular brackets notify an ensemble average equivalent to averaging over many 
periods of the wave. 

The electromagnetic Stokes parameters are defined with respect to two of the observer's 
tetrad vectors e a ,^ and e a (2) as follows [2] 



S 


= Jaf3 


(1) (1) + e (2) e (2) 


(3.29) 


Si 


— Ja/3 


a a 
e (l) e (1) ~ e (2) e (2)_ 


(3.30) 


s 2 


= Ja.fi 


e °(l) e %) + e0l (2) e \l) 


(3.31) 


s 3 


= iJ a/ 3 


ex. (3 ex. (3 
e (l) e (2) _e (2) e (l) 


(3.32) 



It is worth noting that both the Stokes parameters and the components of the polarization 
tensor of an electromagnetic wave can be determined from simple experiments by measuring 
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|£ (1) l 2 + 


|£(2)| 


\c i 2 
|£(i)l - 


1^(2)1 



the intensities of the wave after it passes through devices that transmit only radiation of 
certain polarization [22]. 

Using the definition of the polarization tensor (3.28) we can also rewrite the expressions 
for the Stokes parameters in terms of the electric field components 

12 > , (3.33) 

2 - , (3.34) 

5 2 =< £(i)£(2) + £(i)£(2) > , (3.35) 

5 3 = i < £(i)£(2) - £(i)£(2) > , (3.36) 

where £( n ) = S a e a ^ for n — 1,2. The Stokes parameters can also be expressed in terms of 
the complex functions $ and \1/ since 

l a u 

£* = — y ($ ma + * ™ q) (3 - 37) 

and therefore 

£(D = -^f (* + *), £(2) = (* + *). (3.38) 

To determine how the polarization properties of light vary along the ray we should specify 
the law of propagation for the tetrad vectors. The vectors of both the null tetrad and the 
tetrad e"^ are postulated to be parallelly transported along the light rays. We specify the 
tetrads at the point of observation and thus specify them at every point along the ray. Thus, 
the equations for propagation of the tetrad vectors read 

^ + r^m^ = 0, (3.39) 
of A 

^ + 1^/^ = 0, (3.40) 

de a 

and the same laws hold for vectors l a and n a of the null tetrad. 

It follows from (3.12), (3.27) and (3.41) that the complex amplitude of the electric field 
propagates along the ray according to the law 

D\£ a + $£ a = 0. (3.42) 

For the complex amplitudes $ and ^ one has 

D X <Z> + $<Z> = , D A tf + #tf = 0. (3.43) 
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Finally, from the equations (3.33) - (3.36) it follows that the Stokes parameters propagate 
according to the law 

D\S a + 2dS a = , (3.44) 

where a — 0, 1,2,3. We would like to note that despite suggestive notation the Stokes pa- 
rameters do not form a four-vector, because they do not behave as a vector under coordinate 
system transformations. 

Any stationary or time-dependent axisymmetric gravitational field in general causes a 
relativistic effect of rotation of the polarization plane of an electromagnetic wave [3]. This 
effect was first discussed by Skrotskii [17] and afterwards by many researchers (see, for 
example, [18] and references therein). Recently the effect was studied in [3] where authors 
derived an expression for the angle of rotation of the polarization plane of an electromagnetic 
wave propagating through a weak gravitational field described by metric perturbations h af3 . 
Here we give this expression without derivation which can be found in [3] or in [6]. If is 
the characterizing the orientation of the polarization ellipse, then one has 

—j- = —k k? ' £jpqdqh a p , (3.45) 

where the hat over the spatial indices denotes the projection onto the plane orthogonal to the 
propagation of light ray, for instance, A 1 = P)A> . In the following section we shall integrate 
this equation analytically to obtain the angle of the rotation of the polarization plane for an 
electromagnetic signal propagating in the gravitational field of an isolated system. 

4. SOLUTION TO THE EQUATIONS OF PROPAGATION OF 
ELECTROMAGNETIC SIGNALS 

4.1. Method of Analytical Integration of the Equations of Propagation of 

Electromagnetic Signals 

In this section we describe the method of analytical integration of equations of light prop- 
agation in the field of an isolated gravitating system developed in the series of publications 
[15], [1] and [6]. In subsequent sections this method will be applied to integrate the equa- 
tions of geodesies (3.17) and the equation for Skrotskii effect (3.45) in the space-time with 
the metric given by (2.7) around an isolated gravitating system. 
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We introduce new variables r and as follows 



(4-1) 



where P* is the operator of projection onto the plane perpendicular to k l . If one considers 
unperturbed trajectory of a photon x l N = k % {t — t ) + xo, it is easy to see that the variable 
r characterizing the position of the photon is proportional to time 



T 



N 



t-t*, 



(4.2) 



where t* = kiX l — 1 is the time of the closest approach of the electromagnetic signal to the 
origin of the coordinate system. Since t* is a constant for a particular ray, one has dr = dt. 
This allows to change to the variable r when calculating integrals along the unperturbed 
ray. Vector is the vector from the origin of the coordinate system to the point of the 
closest approach. In terms of the new variables the unperturbed trajectory can be written 
as 

xM=k i r + ?. (4.3) 

Since the vectors k l and £ l are orthogonal to each other, the distance from the point on 
unperturbed trajectory with coordinates r and to the origin of the coordinate system can 
be expressed as 

r N = V^Td* , (4.4) 

where d — |£| is the impact parameter of the unperturbed light-ray trajectory with respect 
to the origin of the coordinate system. 

We introduce the operators of differentiation with respect to r and £ l 



Q. = pi — 
dt ~ ^ W ■ 



(4.5) 



Then for any smooth function F(t, x l ) the following relationships are valid between the 
derivatives in the old and new variables taken on the unperturbed trajectory 



- x=xo+k(t—to) 



= ±F(t*+T,t + kT). 



(4.6) 
(4.7) 



x=xo+k(t— to) 

In the left-hand sides of the equations (4.6) and (4.7) one has to first calculate the derivatives 
and only after that substitute the unperturbed trajectory x = x + k(t — t ) while in the 
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right hand sides one has to first substitute the unperturbed trajectory parameterized by the 
variables r and £ l and then differentiate. Using the equation (4.7) in (4.6) one can obtain 
the expression for spacial derivatives 

~dF(t,x) 

x=xo+k(t—to) 

Equations (4.7) and (4.8) can be used for changing over to the variables r and in the 
equations of light propagation. A useful property of the variables r and £ l is that when one 
calculates integrals along the light rays the following formulae are valid 



{w +k i- k *l) F ^ +T ^ +k ^- < 4 - 8 > 



J §pF(r,t) dr =F(r,t) + C(Z), 



(4.9) 



/ ^nr,t)dT=-^jF(T,t)dT, (4.10) 

where C(£) is a function of £. Equation (4.9) shows that the terms represented as partial 
derivatives with respect to r can be immediately integrated. Equation (4.10) states that one 
can change the order of integration and differentiation which crucially simplifies the problem 
of integration of the equations as we show below. 

Let us consider the geodesic equations (3.17). All terms in the right-hand sides of the 
equations are proportional to the first-order derivatives of the metric perturbations h a p with 
respect to time and spacial coordinates. The metric perturbations, given by the equations 
(2.7), consist of the canonical and gauge- dependent parts. The canonical part (Eqns. (2.2)- 
(2.5)) is expressed as a linear combination of derivatives of different orders with respect 
to spacial variables of functions [F(t — r)/r] (where F(t — r) denotes a mass- or spin-type 
multipole moment of the system). If one changes to the variables r and £ l in the geodesic 
equations (3.17), integrates and changes the order of differentiation and integration in the 
right-hand sides of the equations, then the only two types of integrals that will appear 
in the solutions (for a moment, leave apart the terms produced by the gauge-dependent 
components in the metric perturbations) will be the integrals of the type [F(t — r)/r]'~ 1 l 
and [F(t - r)/r]^ 2 l The gauge-dependent terms, after changing the variables to r and £*, 
will appear in the equations of geodesies under the second-order derivative with respect to 
t and thus can be immediately integrated (cf. Eq. (4.10)). A similar consideration of the 
integration procedure can be done for the equations describing the Skrotskii effect. 
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As it follows from the consideration above, the problem of integration of the equations of 
light propagation reduces to evaluation of integrals [F(t — r)/r]H and [F(t — r)/r]^ where 
F(t — r) denotes multipole moments of different type and order of the gravitating system. 

Let us consider first the contributions to the relativistic effects in propagation of light 
due to the mass-monopole and spin-dipole terms. We neglect the loss of energy and angular 
momentum due to gravitational radiation so that the total mass and angular momentum 
of the system are considered to be conserved. Then it can be shown that the contributions 



under consideration are expressed in terms of integrals [1/r] 



4(l/r) 



-2 



and 



4a(l/r) 



-2] 



4(l/r) 



4a(l/r) 



which can be evaluated. 
The relativistic corrections to the light ray path and the relativistic effect of rotation of 
the polarization plane due to higher order multipole moments (starting with the mass- and 
spin-quadrupole) are expressed in terms of integrals [F(t — r)/r]' _1 ' and [F(t — r)/r]'~ 2 ] along 
the unperturbed light ray, where F(t — r) denotes a multipole moment. To evaluate these 
integrals we introduce the new variable [1] 



V 



= s — t* = t — r(r) = t — V d 2 + t 2 



Then the following relationships hold 
y 2 -d 2 



dr = — ay. 



(4.11) 



(4.12) 



2y 2 y 2 y 

Making use of the new variable y and the equations (4.12) we can express the integrals under 
discussion as follows 



F(t - r) 



F(t - r) 



hi] 



y 

-I 



F(t* + {) 



c 



y v 




F(t* + C) 

c 



d 2 



y v 

i r Fit* + o 



— oo — oo 



Tj 

— oo — oo 



c 



d(drj , 



(4.13) 



(4.14) 



where (, rj are new variables of integration, and t* is the time of the closest approach of 
photon to the origin of the coordinate system. It is worth to note that the time t* depends 
on the choice of coordinate system and therefore in general has no physical meaning. 

An important property of the integrals (4.13) and (4.14), expressed in terms of the variable 
y, is that they depend on r and £ l only through the upper limits of integration (since 
y = t — \J d 2 + t 2 ) and the square of the impact parameter d 2 in the prefactor of the second 
integral in (4.14). 
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Differentiating (4.13) with respect to £ l and r yields 





~F(t-r)~ 


'""! 




r 




[ 


~F(t - r)~ 


'-"1 




r 





F(r + y) $y = -Fit* + y) 4 ln(-y) = ^F(t- r) 

y yr 



F(t* + y) 



dry 



F(t* + y) 



F(t-r) 



(4.15) 



(4.16) 



y y 

The last result also follows from the equation (4.10). Thus derivatives of the integral [F(t — 
r)/r]H with respect to either or r can be expressed in terms of the integrand. 
Differentiating the (4.14) with respect to £ l multiple times yields 




F(t* + C) 

c 



d(drj 



(4.17) 
(4.18) 



da 



ijk 



F(t 






r) 




r 




—00 















F(t* + C) 

c 

F(t- 



d(drj 



(4.19) 



In the last expression all terms contain at least a first order derivative of the integral [F(t — 
r)/r]l _1 l with respect to £ l . Above it was shown (Eq. (4.15)) that in this case the integration 
is eliminated. Therefore the third and higher order derivatives of [F(t — r)/r] '~ 2 1 with respect 
to £ l can be expressed in terms of F(t — r) and one does not have to perform integration. 
Explicitly evaluating the derivatives in (4.19) using (4.15) one gets 



d, 



ijk 



F(t-r) 



pij^k _|_ pjk^i _|_ pik^j 



y 2 r 



F(t 



(4.20) 



+ 



egg 



\y rj 



F(t - r) - F(t - r) 



In the equations for trajectory of a photon and gravity-induced rotation of the polarization 
plane the integrals [F(t — r)/r]H and [F(t — r)/r]^ 2 ^ are differentiated with respect to r 
and e and it can be shown that in all terms differentiation eliminates integration. In the 
subsequent sections we write out solutions to the equations formally in terms of integrals 
[F(t — r)/r]H and [F(t — r)/r}^~ 2 \ but one should bear in mind that performing differen- 
tiations using the formulae (4.9), (4.15) and (4.19) the integrals can be expressed in terms 
of integrands. 
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■ rr {w i ~ Vw ) , (4.21) 



4.2. Solution to the Equations for Trajectory of a Photon 

Using the formulae (4.7) and (4.8) and rewriting the equation (3.17) in new variables one 
obtains 

dr 2 

where all functions in the right-hand side are taken on the unperturbed light-ray trajectory 
(before performing differentiation). The gauge functions w a have not yet been specified 
which means that equations (4.21) are gauge-invariant. Substituting into (4.21) the met- 
ric perturbations expressed in terms of r and (Eqns. (A.1)-(A.9) in Appendix A) and 
integrating once with respect to r one obtains [6] 

x i (T) = k* + Z\T,H), (4.22) 

where 

S*(r, = E l (r, + (r, £) + (r, £) , (4.23) 

(G) (M) (5) 



S< (r, £) =<9 r [(^ - fcV°) + K - Vw )] , 



S* (r, £) =2M 



(M) 



s. 1 - 

r 

oo i p 



=2 p=0 g=0 



r 

Z! 



(4.24) 
(4.25) 



C l (l-p,p-q,q)H(2-q)x 



1 - 



P - g 

z 



1 - 



Z- 1 



k< ai ... a „d t 



<ai ...a p L 'ap-)-i . ..a^ >^7- 



Xt q \t-r) 



;=2 p=0 



1 ^ ^ J ZCj< ai ... ap 9 ap+1 ... ai > 



2p 



k< ai ...a p -ida p ...a l _ 1 > 
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E l (r, £) = 2kj€ jba S b 



(S) 



-2a 



tibaSb 



(4.26) 



oo i— 1 p 



4*A E E E VriTr ^- 1 (/ " p ~ x ' p " g ' g) F(2 ~ q) x 

Z=2 p=0 g=0 ^ >' 



1 - 



P-9 



£ / ■-<ai...a P L/ a p +i...a;_i> ty r 



' k<ai...a v d( 

4 (4 - fc a 4*) e E TTTryr^- 1 (' - p - * 



+ 



Z=2 p=0 



P \ 

I _Y / <°i— °p °p+i— °J-i> 



^^(i-r) 



+ 



oo i— 1 



Z=2 p=0 



(( + 1)! 



Here (r, £) represents the gauge-dependent perturbations, S 8 (r, £) and S 8 (r, £) are the 

(G) ' ' (M) (S) 

perturbations due to the mass and spin multipole moments, correspondingly. The Heaviside 
function H(p — q) is defined by the expression (1.18) and Ci(l—p,p — q, q) are the polynomial 
coefficients (1.17). The gauge functions (p a were introduced as follows: we collected in the 
equations (4.21) all terms with second and higher order derivatives with respect to r and 
equated them to k l p° — p l . By introducing the functions (p a we single out the terms that can 
be immediately integrated. We do not separate k l p Q — tp l into <p° and ip 1 , such separation 
is not unique while the functions k l p° — if 1 are uniquely defined (Eqns. (B.4)-(B.6) in the 
Appendix B). 

The gauge functions w a can be chosen arbitrarily. For the reasons discussed in Section 
2 we choose w a which make our coordinate system ADM-harmonic, that is satisfying both 
ADM and harmonic gauge conditions. These functions are given by the expressions (B.2) 
and (B.3). 

The second integration of (4.21) yields the expressions for trajectory of the photon 



x t (r)=x i N + AE l +AE t +AS ! , 

(G) (M) (S) 



(4.27) 
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where 



(G) (G) (G) 
(M) (M) (M) 

AS ! eS ! (r,0-S J (r ,€) • 

(S) (S) (S) 



(4.28) 



The term 



(G) 



(4.29) 



represents the gauge- dependent part of the trajectory's perturbation, and the physically 
meaningful perturbations due to the mass and spin multipoles are given by 



~ 1 (r, £) =2M 



(M) 



[-2] 



2.MF 



+ 



oo I p 



£=2 p=0 q=0 



Z! 



_ p-g \ / , p-q , , 



2$-*>(t-r) 



oo Z— 1 
Z=2 p=0 



^ ]_ / ki<ai...a p d a , p+1 ...a l > 



(4.30) 
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~ l (r, £) =2kj€ jba S b 



(S) 



"1 [-2] 



r 

oo I— 1 p 



{-i) l+ p-n 



4i 

r 



(4.31) 



4*A £ £ £ (/ • n! ,( - < ; - /' 1 • /' - «■ '/ > H ( 2 - </) x 



Z=2 p=0 q=0 

P-9 



I _ ^ ) K <ai...a p O ap+1 ...a l _ 1 : 



+ 



Z=2 p=0 



P 



I _ I I ^<°i — a p 9 ap+1 ...ai- 1 > 



e iba S^Jt - r) 



+ 



00 i_1 r_iV+p/ 

^££tzw <7, - i(z " p " 1 ' p)x 



(Z + l)! 

,^2 p=0 V ; 
k<a 1 ...a p da p+ i...ai- 1 > 



4.3. Solution to the Equations for Rotation of the Polarization Plane 

We rewrite the equation (3.45) in terms of the variables r and £ l using the relation- 
ship (4.8) and substitute the metric (2.7) expressed in terms of new variables (Eqns. (A.l)- 
(A.9)): 

d f r = \k a Ve m d q Kp- + \k j e m d qT w p ■ (4.32) 
Integrating with respect to r and changing the order of differentiation and integration yields 

= 0(G) + 0(M) + 0(S) + 00 , (4.33) 

where 0o is the constant angle characterizing the initial orientation of the polarization ellipse 
in the plane formed by the vectors em and em. The terms (G) , (M) and (5) describe 
the contributions due to gauge-dependent terms, terms depending on mass- and spin-type 
moments correspondingly. 

The gauge-dependent part of the Skrotskii effect can be obtained immediately, since the 
gauge dependent terms appear in the equation (4.32) under the derivative with respect to 
t (cf. Eq. (4.9)). The integration yields 

1. A 



0(G) = ^Wejpqdq (w p + x p ) , 



(4.34) 
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where the gauge functions w l and x % are given by the equations (2.15) and (B.7). The 
gauge functions \ % were introduced by collecting in the equation (4.32) all terms that can 
be eliminated by a gauge transformation. 

The Skrotskii effect due to the mass and spin multipoles of the isolated system is given 

by 

(4.35) 



1=2 p=0 



/! 



I- 1 



k 8 

fb <ai...a p L/ a p+ i...a;>j 



and 



< ^~ > ^ /" 1 \ l-\-T) / 

feW= 2 EEVTTT c ' (; - p ' p »( 1 -7) 



(4.36) 



l + iJ(/-l) 1- 



2p 



rv <ai...a p '-'ap + i...a l > 



[-i] 



Just like in the case with the equations for light-ray path (4.25)-(4.26) and (4.30)-(4.31) 
differentiation eliminates integrals in the right-hand sides of the equations (4.35)-(4.36) as 
it was shown in Section 44.1. 

5. OBSERVABLE RELATIVISTIC EFFECTS 

In this section we give the expressions for the observable relativistic effects of time delay 
and gravitational deflection of light. The expressions for another observable effect - rotation 
of the polarization plane were essentially given in Section 44.3. 

5.1. Time Delay 

In the paper [1] the following expression was obtained for the time of propagation of the 
electromagnetic signals through the gravitational field of an isolated source from the emitter 
to the observer 



t - t =\x - x \ + A(r,r ) , 
A(r, r ) = A (r, r ) + A (r, r ) + A (r, r ) , 

(G) (M) (S) 



(5.1) 
(5.2) 
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where Xq = (t , x Q ) and x a = (t, x) are the coordinates of the points of emission and 
observation of the signal, A (r, r ), A (r, r ) and A (r, r ) are functions describing the delay 

(G) (M) (S) 

of the electromagnetic signal due to the gauge-dependent terms, mass and spin multipoles of 
the gravitational field of the isolated system correspondingly These functions are expressed 
in terms of the perturbations of the signal's trajectory as 



A (r, r ) = 

(G) 

A (r,r ) = 

(M) 



A ; 



A (r, r ) = - h 

(S) 



S< (r,£)-H* (to,*) 

(G) (G) 
(M) (M) 

S< (r,0-2* (t-0,0 

(S) (S) 



(5.3) 
(5.4) 
(5.5) 



We note that the expressions 5.3 give the effect with respect to the coordinate time which 
has to be converted to the time measured by the observer. It is shown in [1] that if the 
observer's velocity is negligible with respect to the global ADM-harmonic coordinate system 
the relationship between the coordinate time and proper time T of the observer is given by 

) , (5-6) 



T = 



where t\ is the initial epoch of observation. In most cases the distance r is much grater than 
M. and coordinate time coincide with the proper time of the observer. If observer is moving 
with respect to the global coordinate system, an additional Lorentz transformation to the 
rest frame of the observer has to be performed. This case was considered in [6]. 

5.2. Deflection of Light 



The expression for observable vector towards the source of light calculated with respect 
to the observer's reference frame was derived in [1] and reads 

s l (r, £)=K* + a 1 (r, £) + (r, £) - ? (r , £) + Y (r, £) , (5.7) 

where 



(5.8) 



|x - x | 

is the unit vector in the direction "observer - source of light" , defined in the global coordinate 
system; 



a 1 (T,£) = -PjE> 



(5.9) 
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is the vector describing the angle of light deflection; 



P (r, 





3* (r, |) + 3* (r, £) 

_(M) (S) 




x - x 





(5.10) 



is relativistic corrections, introduced by the relationship 

fc* = -/r-/P(T,$)+/F(T ,$)- (5.11) 

The term 



y(r,0=-^'^ T (t,x) (5.12) 



appears as a result of transformation from the the global ADM-harmonic system to the local 
frame of the observer and describes the perturbations of the observer's coordinate system 
due to the gravitational waves emitted by the isolated system. It was assumed that the 
observer is at rest with respect to the global ADM-harmonic system. In the case when 
the observer is moving with respect to the global system one has to perform an additional 
Lorentz transformation to the rest frame of the observer. 

APPENDIX A: THE METRIC TENSOR IN TERMS OF VARIABLES £ AND r 

The expressions for perturbations of the metric tensor (2.2) - (2.5) in variables £ and r 
have the form as follows 

9 KA oo I p 

1=2 p=0 q=0 



1=2 p=0 q=0 



(S) (M) 



(A.2) 



= W- + C 1 ' > ( A - 3 ) 



oo 1—2 p oo /— 2 p 



Z=2 p=0 g=0 (M) Z=2 p=0 g=0 (S) 
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where 
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All quantities in the right side of expressions (A.5)-(A.9), which are explicitly shown as 
functions of x l , r = \x\ and t, must be understood as taken on the unperturbed light-ray 
trajectory and expressed in terms of £ l , d — r and t* in accordance with the equations 
(4.2), (4.4). For example, the ratio X^~ q \t — r)/r in equation (A.5) must be understood as 



2<P-«)(t _ r ) l%- q) (t* + r- V^Td?) 



-(p-q) i 
- -A 



r ' A 2 + d 2 

and the same replacement rule is applied to the other equations. 



(A.10) 



APPENDIX B: GAUGE FUNCTIONS 

Gauge functions w a , generating the coordinate transformation from the canonical har- 
monic coordinate system to the ADM-harmonic, are given by equations (2.14), (2.15). They 
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transform the metric tensor as follows 



(B.l) 



where h^S" is the canonical form of the metric tensor in harmonic coordinates given by 
equations (2.2)-(2.5) and h a p is the metric tensor given in the ADM-harmonic coordinates 
by equations (2.16)-(2.19). 

The gauge functions taken on the light-ray trajectory and expressed in terms of the 
variables and r can be written down in the form 
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where 1iqq 9 (-u, r, £), h^ 9 (u, t, £) and h^ 9 (tt, r, £) are defined by the Eqns. (A. 5), (A. 6) and 

(M) ' (M) (S) 

(A. 7) after making use of the substitution i* — > it. 

Linear combination — of the gauge-dependent functions tp a introduced in equation 
(4.24) is given by the expressions 



(B.4) 
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Gauge-dependent term generated by equation (4.34) for the rotation of the plane of polar- 
ization of electromagnetic wave is a pure spatial vector x % that can be decomposed in two 
linear parts corresponding to the mass and spin multipoles: 



"X. X(m) "X- 
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where 
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